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Optical  absorption  in  molecular  crystals 
from  time-dependent  density  functional  theory 

Final  report  for  the  European  Office  of  Aerospace  Research  &  Development 

Grant  No.  FA9550-15-1-0290 

Leeor  Kronik1 

1  Department  of  Materials  and  Interfaces,  Weizmann  Institute  of  Science,  Rehovoth  76100,  Israel 

We  address  a  major  challenge  for  computational  materials  science  based  on  density  functional 
theory,  by  showing  that  fundamental  gaps  and  optical  spectra  of  molecular  solids  can  be  predicted 
quantitatively  and  non-empirically  within  the  framework  of  time-dependent  density  functional  the¬ 
ory  (TDDFT),  using  the  recently-developed  optimally-tuned  screened  range-separated  hybrid  (OT- 
SRSH)  approach.  In  this  scheme,  the  electronic  structure  of  the  gas-phase  molecule  is  determined  by 
optimal  tuning  of  the  range-separation  parameter  in  a  range-separated  hybrid  functional.  Screening 
and  polarization  in  the  solid-state  are  taken  into  account  by  adding  long-range  dielectric  screening  to 
the  functional  form,  with  the  modified  functional  used  to  perform  self-consistent  periodic-boundary 
calculations  for  the  crystalline  solid.  We  provide  a  comprehensive  benchmark  for  the  accuracy  of 
our  approach  by  considering  the  X23  set  of  molecular  solids  and  comparing  results  obtained  from 
TDDFT  with  those  obtained  from  many-body  perturbation  theory  in  the  GW-BSE  approximation. 

We  additionally  compare  results  obtained  from  dielectric  screening  computed  within  the  random- 
phase  approximation  to  those  obtained  from  the  computationally  easier  many-body  dispersion  ap¬ 
proach  and  find  that  this  influences  the  fundamental  gap  but  there  is  little  effect  on  the  optical 
spectra.  We  therefore  believe  that  the  method  is  robust  and  can  be  used  for  studies  of  molecular 
solids  that  are  typically  outside  the  reach  of  computationally  more  intensive  methods.  The  work  is 
now  being  written  up  for  publication. 


I.  INTRODUCTION 

The  electronic  and  optical  properties  of  molecular 
solids  have  recently  attracted  significant  attention,  pri¬ 
marily  in  the  context  of  optoelectronic  devices  based 
on  small  molecules  (see,  e.g.,  Refs.  1-4).  In  particular, 
there  is  on-going  interest  in  identifying  small-gap  organic 
molecules  for  high-performance,  low-cost,  or  enhanced- 
stability  optoelectronic  devices  based  on  solids  synthe¬ 
sized  from  these  molecules  (see,  e.g.,  Refs.  5-9,  for  some 
recent  overviews) .  Theory  can  and  should  play  an  impor¬ 
tant  role  in  such  investigations,  as  it  can  clarify  the  prop¬ 
erties  of  existing  molecular  solids  and  point  out  promising 
new  ones  [10-13]. 

Electronic  properties,  such  as  the  band  structure  in 
general  and  the  transport  gap  in  particular,  and  opti¬ 
cal  properties,  such  as  optical  absorption  in  general  and 
the  optical  gap  in  particular,  are  excited-state  proper¬ 
ties.  For  inorganic  solids,  these  properties  have  long 
been  calculated  using  many-body  perturbation  theory 
(MBPT). [14-16]  In  MBPT,  Dyson’s  equation  is  often 
solved  using  Hedin’s  GW  approximation, [17]  where  G  is 
the  one-particle  Green  function  and  W  is  the  dynami¬ 
cally  screened  Coulomb  interaction.  [14,  18]  The  Bethe- 
Salpeter  equation  (BSE)  for  the  two-particle  Green  func¬ 
tion  is  then  solved  approximately  to  predict  optical 
properties.  [14,  19,  20]  In  recent  years,  the  GW-BSE  ap¬ 
proach  has  been  increasingly  applied  to  molecular  solids, 
yielding  many  important  insights  (see  Ref.  13  for  a  re¬ 
cent  overview).  Unfortunately,  such  GW-BSE  calcula¬ 
tions  can  be  quite  expensive  and  complicated,  limiting 


our  ability  to  use  them  routinely,  especially  in  the  con¬ 
text  of  high-throughput  calculations  for  new  materials. 

Density  functional  theory  (DFT),  in  both  its  time- 
independent  [21,  22]  and  time-dependent  [23-26]  forms, 
suitable  for  ground  and  excited  state  properties,  respec¬ 
tively,  is  much  more  efficient  computationally.  How¬ 
ever,  common  approximations  to  time-dependent  DFT 
(TDDFT)  are  known  to  fail  in  the  solid-state  limit 
[14,  27].  For  molecular  solids  in  particular,  key  quan¬ 
tities,  such  as  the  transport  gap,  the  optical  gap,  and  the 
exciton  binding  energy  (i.e.,  the  difference  between  the 
two  gaps),  are  often  in  qualitative  or  gross  quantitative 
error.  [13] 

Recently,  Refaely- Abramson  et  al.  have  suggested  the 
optimally-tuned  screened  range-separated  hybrid  (OT- 
SRSH)  functional  as  a  means  for  quantitative  DFT- 
based  prediction  of  excited-state  properties  in  molecu¬ 
lar  solids.  [28,  29]  In  this  approach,  one  first  computes 
the  underlying  gas-phase  molecule  using  an  asymptot¬ 
ically  correct  range-separated  hybrid  (RSH)  functional, 
in  which  an  optimal  range-separation  parameter  is  de¬ 
termined  non-empirically,  [30-33]  based  on  satisfaction  of 
the  ionization  potential  theorem.  One  then  uses  the 
same  range-separation  parameter  in  the  solid-state  en¬ 
vironment,  while  accounting  explicitly  for  solid-state  po¬ 
larization  by  screening  the  asymptotic  potential  with  a 
non-empirical  dielectric  constant.  [28,  29] 

While  preliminary  results  obtained  with  the  OT-SRSH 
method  have  shown  excellent  agreement  with  GW-BSE 
data,  two  important  questions  remain.  First,  results  have 
been  reported  to-date  only  for  only  two  molecular  solids  - 
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pentacene[28,  29]  and  quinacridone,[34]  an  air- st able  pen- 
tacene  derivative,  with  the  optical  absorption  spectrum 
spectra  computed  only  for  the  former.  The  validity  of 
the  OT-SRSH  approach  across  a  wider  range  of  molecu¬ 
lar  crystals  is  therefore  in  need  of  demonstration.  Second, 
previous  OT-SRSH  calculations  have  used  the  dielectric 
constant  obtained  within  the  random-phase  approxima¬ 
tion  (RPA)  for  facilitating  comparison  to  MBPT  data. 
However,  this  step  can  itself  be  expensive  and  it  remains 
to  be  seen  whether  sufficiently  accurate  results  can  be 
obtained  from  more  simple  methods  for  determining  the 
dielectric  constant. 

In  this  article,  we  address  both  questions  by  assess¬ 
ing  the  accuracy  of  the  OT-SRSH  approach  for  trans¬ 
port  gaps  and  optical  absorption  spectra  across  the  X23 
set  of  molecular  solids.  [35,  36]  This  set  comprises  crystals 
based  on  small-  to  medium-sized  organic  molecules,  pos¬ 
sessing  a  variety  of  weak  inter-molecular  interactions  and 
different  degrees  of  solid-state  polarization.  It  therefore 
provides  a  strict  benchmark  for  OT-SRSH  capabilities. 
We  further  compare  results  obtained  using  an  RPA-based 
dielectric  constant  with  those  obtained  using  many-body 
dispersion  (MBD).  Within  the  RPA,  we  find  our  OT- 
SRSH  results  to  be  in  very  good  agreement  with  those 
obtained  from  GW  for  quasi-particle  gaps  and  from  GW- 
BSE  for  the  optical  spectrum.  We  further  find  that  using 
MBD-based  dielectric  screening  results  in  larger  devia¬ 
tions  for  quasi-particle  gaps  but  has  an  essentially  neg¬ 
ligible  effect  on  the  optical  absorption,  allowing  for  an 
inexpensive  yet  non-empirical  prediction  of  optical  prop¬ 
erties. 


II.  THEORETICAL  AND  COMPUTATIONAL 
APPROACH 

A.  Optimally-tuned  range-separated  hybrid 
functionals 

In  the  range-separated  hybrid  (RSH)  method,  the 
Coulomb  interaction  is  range-split.  Here,  we  use  the 
range-separation  scheme  suggested  by  Yanai  et  al, [37] 
which  is  based  on  the  identity: 

1  a  +  /3erf(7r)  ^  1  —  [a  +  f3ei-i("/r)\ 
r  r  r 

where  r  is  the  inter-electron  coordinate  and  a,  /?,  7  are 
parameters.  The  full  1/r  repulsion  is  used  for  the  Hartree 
and  correlation  terms,  but  the  two  terms  on  the  right- 
hand  side  of  Eq.  (1)  are  treated  differently  in  the  compu¬ 
tation  of  the  exchange  term.  The  first  term  is  treated  us¬ 
ing  exact  (i.e.,  Fock)  exchange,  whereas  the  second  term 
is  treated  using  local  or  semi-local  exchange.  This  leads 
to  the  following  expression  for  the  exchange-correlation 
energy,  Exc: 

£*sh  =  aE XSXR  +  (1  -  a)E gRAx  +  (a  +  /1)FXLR 

+  (1  —  a  —  /3)Ed  fax  +  -E’dfAc;  (2) 


where  the  super-scripts  ‘SR’  and  ‘LR’  denote  short-range 
and  long-range  contributions,  respectively,  and  the  sub¬ 
scripts  ‘xx’,  ‘DFAx’  and  ‘DFAc’  denote  Fock-like  ex¬ 
act  exchange,  approximate  (semi-)local  exchange,  and 
approximate  (semi-)local  correlation,  respectively.  [38] 
Equation  (2)  reveals  that  the  parameter  a  dictates  the 
amount  of  Fock-like  exchange  in  the  short  range  (r  — >  0) 
and  the  parameter  sum  a  +  (3  determines  the  amount  of 
Fock-like  exchange  in  the  long  range  (r  — »  00).  The  two 
limits  are  smoothly  interpolated  using  the  error  function, 
with  7  being  the  range-separation  parameter,  i.e.,  I/7 
corresponds  to  a  typical  length  denoting  the  transition 
from  SR  to  LR. 

In  order  to  turn  Eq.  (2)  into  a  practical  functional,  one 
needs  to  choose  the  approximate  (semi-)local  exchange- 
correlation  functional  and  set  the  parameters  a,  /3,  and 
7.  To  proceed  without  introducing  empiricism,  typical 
choices  for  the  (semi)-local  functional  would  be  the  local 
density  approximation,  LDA,[39]  or  the  Perdew-Burke- 
Ernzerhof  (PBE)[40]  form  of  the  generalized-gradient  ap¬ 
proximation  (GGA).  In  some  cases,  the  fraction  of  SR 
Fock  exchange,  a,  can  be  determined  from  first-principles 
based  on  satisfaction  of  piecewise  linearity  in  fractional 
DFT  [38,  41],  but  this  is  not  always  possible.  [42]  For  a 
wide  variety  of  organic  molecules  a  universal  value  of  0.2 
has  been  found  to  be  useful  (see,  e.g.,  Refs.  28,  38,  41- 
43).  This  value  is  used  in  this  work  throughout.  For  any 
choice  of  a,  the  condition  a  +  f3  =  1  guarantess  100%  of 
LR  Fock  exchange  and  therefore  the  correct  asymptotic 
potential  in  the  gas  phase.  [30,  38,  41]. 

In  many  popular  RSH  functionals,  the  range- 
separation  parameter,  7,  is  given  a  universal  value  based 
on  fitting  against  an  appropriate  data  set.  [31,  37,  43- 
45]  In  the  optimal  tuning  (OT)  scheme,  7  is  system- 
dependent,  but  still  chosen  non-empirically.  For  gas- 
phase  systems,  it  is  obtained  by  satisfying  the  ioniza¬ 
tion  potential  (IP)  theorem, [46-49]  which  states  that  for 
the  exact  functional  the  energy  of  the  highest  occupied 
molecular  orbital  (HOMO)  is  equal  and  opposite  to  the 
ionization  potential,  i.e.,  IP  =  — chomo ■  Often,  this 
condition  is  demanded  simultaneously  for  the  system  in 
both  its  neutral  and  anionic  state  (where  the  ionization 
potential  corresponds  to  the  electron  affinity  of  the  neu¬ 
tral).  [30,  31,  50]  For  any  choice  of  a,  the  optimal  tuning 
scheme  then  involves  the  minimization  of  a  target  func¬ 
tion,  J(q;a),  defined  by: 

J2(7;a)=  E  [IP^(i)+sY(i)?-  (3) 

i—N,N+l 

Minimizing  this  target  function  has  been  shown  to  be 
equivalent  to  enforcing  piecewise  linearity  [33,  34,  41,  51- 
53],  resulting  in  an  accurate  prediction  of  the  ionization 
potential  and  the  electron  affinity  directly  from  the  en¬ 
ergy  levels  of  the  highest  occupied  and  lowest  unoccupied 
orbitals,  respectively.  [31,  33,  54] 

Refaely- Abramson  et  al.  have  suggested  that  the  OT- 
RSH  scheme  can  be  extended  to  molecular  solids  by  using 
a  screened  range-separated  hybrid  (SRSH).[28]  Briefly,  in 
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this  approach  one  first  selects  a  and  7  for  the  gas-phase 
molecule  as  described  above.  One  then  notes  that  in  the 
gas- phase  the  asymptotic  potential  is  —  1  /r  but  in  the 
solid  state  it  is  —  l/(er),  where  e  is  the  scalar  dielectric 
constant.  Therefore  (3  is  re-adjusted  to  reflect  this  screen¬ 
ing,  by  demanding  that  a+/3  =  1/e  instead  of  a+/3  =  1  as 
in  the  gas-phase.  The  resulting  screened  RSH  functional 
is  then  applied  to  the  molecular  solid.  [13,  28,  29,  34] 

B.  The  X23  set  of  molecular  solids 


functional, [43],  which  is  based  on  Eq.  2  with  a=0.2  and 
PBE  correlation  and  short-range  exchange  components, 
as  implemented  in  the  Q-CHEM  code  (version  4.3)  [55], 
but  with  the  range-separated  parameter  7  optimally- 
tuned  per  system,  rather  than  fixed  to  its  default  value. 
Optimization  proceeded  via  minimization  of  the  target 
function  J  given  in  Eq.  (3),  i.e. ,  both  neutral  and  an¬ 
ion  were  considered,  except  for  molecules  exhibiting  an 
unbound  LUMO,  where  only  the  neutral  form  was  con¬ 
sidered.  The  all-electron  cc-pVTZ  basis  set  was  used 
throughout  for  all  atoms. 


For  evaluating  the  accuracy  of  the  above  approach  in 
a  systematic  manner,  we  consider  quasi-particle  (QP) 
gaps  and  optical  absorption  spectra  for  the  set  of  23 
non-covalently  bound  molecular  solids,  known  as  the 
X23  set. [35]  A  schematic  diagram  displaying  all  molec¬ 
ular  entities  considered  in  this  set  is  given  in  Fig¬ 
ure  1.  The  molecules  used  are  small-  to  medium-sized 
organic  molecules  that  can  be  grouped  into  four  subsets 
based  on  their  chemical  identity:  Open-cyclic  aliphatic 
molecules  (carbon  dioxide,  ammonia,  acetic  acid,  suc¬ 
cinic  acid,  cyanamide,  ethyl  carbamate,  oxalic  acid  (in 
both  a  and  /3  polymorphs),  urea,  and  formamide),  cyclic 
aliphatic  molecules  (adamantane,  hexamine,  trioxane, 
and  1,4-cyclohexane-di-one),  cyclic  aromatic  molecules 
(benzene,  naphthalene,  and  anthracene),  and  hetero¬ 
cyclic  aromatic  molecules  (cytosine,  uracil,  triazine,  imi¬ 
dazole,  pyrazine,  and  pyrazole).  In  crystalline  solid  form, 
these  molecules  are  weakly  bound,  typically  through  El- 
bonding,  7r  —  7r  stacking,  van  der  Waals  interactions,  etc. 
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FIG.  1:  (Color  online)  Chemical  structures  of  all  organic 
molecules  present  in  the  X23  molecular  crystal  set. 


C.  Computational  Details 

All  gas-phase  OT-RSH  calculations  presented  in 
this  work  were  based  on  the  LRC-wPBEO  RSH 


All  solid-state  OT-SRSH  calculations  were  carried 
out  using  PARATEC  (revision  499), [56]  modified  a 
pseudopotential-planewave  code.  Here,  short-range  LDA 
exchange  was  used,  together  with  LDA  correlation,  again 
with  a=0.2  throughout.  Differences  in  tuning  based 
on  LDA  or  PBE  were  found  to  be  insignificant.  See 
Refs.  28,  29  for  more  implementation  details.  LDA-based 
Troullier-Martins  [57]  norm-conserving  pseudopotentials, 
adapted  from  the  ABINIT  website,  [58]  were  used  for  all 
atoms.  An  energy  cutoff  of  816  eV  was  used  throughout. 
Two  different  methods  were  used  to  evaluate  the  scalar 
dielectric  constant  needed  for  the  determination  of  /3  in 
the  solid-state  calculations.  In  one,  we  used  the  random 
phase  approximation  (RPA),  as  used  in  the  Go  Wo  cal¬ 
culations  elaborated  below.  In  the  other,  we  used  the 
many-body  dispersion  (MBD)  approach.  [59] 

For  comparison  purposes,  all  molecular  solids  were 
also  computed  using  a  standard  one-shot  perturbative 
Go  Wo  calculation,  [13,  18]  based  on  the  DFT  eigenval¬ 
ues  and  eigenvectors  obtained  from  am  LDA  calculation 
within  PARATEC.  We  used  a  generalized  plasmon  pole 
model  [18],  implemented  within  the  BerkeleyGW  pack¬ 
age  (trunk  version,  revision  6539). [60]  This  approach  has 
previously  been  established  as  a  quantitatively  useful  tool 
for  the  study  of  molecular  solids  (see  Ref.  13,  and  refer¬ 
ences  therein). 

Optical  spectra  in  the  solid  state  were  computed  using 
TDDFT  with  the  LDA  and  the  OT-SRSH  functional,  as 
well  as  with  the  Bethe-Salpeter  equation  (BSE)  based 
on  the  G0W0  output.  Both  TDDFT  and  BSE  calcula¬ 
tions  were  performed  using  the  BerkeleyGW  package  [60] , 
modified  to  include  TDDFT,  with  incident  light  polar¬ 
ization  averaged  over  the  main  unit-cell  axes.  The  kernel 
was  calculated  on  a  coarse  wavefunction  grid,  then  in¬ 
terpolated  to  a  fine  grid  using  the  interpolation  scheme 
suggested  by  Rohlfing  and  Louie.  [61]  We  used  a  slightly 
shifted  grid  to  generate  the  transition  matrix  elements  in 
the  dielectric  function,  using  a  velocity  operator  to  ap¬ 
proximate  an  incident  light  along  a  specific  direction.  [61] 
Grid  shift  directions  were  set  along  the  a,  b ,  and  c  unit¬ 
cell  axes.  For  large  crystals,  we  only  used  one  grid  in  the 
TDDFT  calculations  as  they  are  easier  to  converge  with 
respect  to  k-point  sampling. 
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III.  RESULTS  AND  DISCUSSION 
A.  Fundamental  gap 

We  begin  our  benchmark  evaluation  by  comparing  the 
fundamental  gaps  computed  using  OT-SRSH  and  GqWq 
methods  for  all  solids  in  the  X23  set.  To  prevent  cloud¬ 
ing  the  comparison  by  a  different  treatment  of  the  dielec¬ 
tric  constant,  all  OT-SRSH  results  presented  are  based 
on  the  dielectric  constant  obtained  within  RPA,  except 
in  the  last  part  of  this  section,  where  comparison  with 
MBD-based  results  is  explicitly  made.  With  both  OT- 
SRSH  and  GW,  we  computed  the  fundamental  gap  as  the 
energy  difference  between  the  highest  occupied  and  low¬ 
est  unoccupied  state  at  the  T  point  of  the  Brillouin  zone, 
as  the  inter-molecular  orbital  hybridization  and  therefore 
band  dispersion  are  very  small. 

All  computed  fundamental  gaps  are  given  in  Table  I, 
where  they  are  additionally  compared  to  LDA-computed 
gaps.  It  is  readily  observed  that  the  OT-SRSH  gaps  agree 
very  well  indeed  with  the  GW  ones.  The  deviation  be¬ 
tween  the  gaps  computed  with  both  methods  is  summa¬ 
rized  graphically  in  Fig.  2.  The  differences  are  usually  0.2 
eV  at  most,  with  a  mean  absolute  deviation  of  only  0.15 
eV.  Only  two  solids  (hexamine  and  imidazole)  exhibited 
a  somewhat  larger  deviation  of  0.3  eV  and  only  one  solid 
(pyrazole)  exhibits  a  larger  deviation  of  0.5  eV.  Not  sur¬ 
prisingly,  these  gaps  are  substantially  larger  than  those 
obtained  with  LDA,  which  is  well-known  to  underesti¬ 
mate  fundamental  gaps  in  general.  [14,  18]  A  comparison 
of  the  solid-state,  OT-SRSH-computed  fundamental  gaps 
with  the  gas-phase,  OT-RSH  computed  ones  (not  shown 
for  brevity)  reveals  that,  as  expected,  the  gas-phase  gaps 
are  substantially  larger  than  the  solid-state  ones  (by  as 
little  as  0.8  eV  and  as  much  as  4.8  eV).  This  reflects 
the  well-known  phenomena  of  polarization-induced  gap- 
renormalization  in  molecular  solids,  [62]  which  is  clearly 
captured  in  the  OT-SRSH  scheme  [28]  but  is  known  to 
be  absent  in  standard  functionals  [13,  63]. 

The  above  results  establish  dielectric  screening  as  key 
to  accurate  treatment  of  molecular  solids.  Indeed,  there 
is  growing  recent  interest  in  solid-state  screening  as  an 
ingredient  in  the  construction  of  density  functionals  in 
general  (see,  e.g.,  [29,  64-68]).  We  also  note  that  specifi¬ 
cally  for  molecular  solids,  there  is  growing  interest  in  em¬ 
bedding  a  range-separated  hybrid  molecular  calculation 
within  a  polarizable  continuum  model  to  mimic  solid- 
state  effects  (see,  e.g.,  [69,  70]).  Our  approach  represents 
a  full  solid-state  calculation,  allowing  for  polarization  ef¬ 
fects  while  still  capable  of  capturing  inter-molecular  dis¬ 
persion  effects  [28]. 

B.  Optical  Absorption  Spectra 

We  next  compare  the  optical  absorption  spectra  calcu¬ 
lated  using  TD-OT-SRSH  with  those  obtained  from  the 
GqWo-BSE  methods.  As  we  are  interested  in  allowed 


TABLE  I:  Fundamental  gaps  of  the  X23  set  of  molecular  solids 
(in  eV),  calculated  using  LDA,  OT-SRSH  and  Go  Wo. 


Molecular  Solid 

Fundamental  Gap  (Eg) 

LDA  OT-SRSH  G0IFo 

Carbon  dioxide 

6.4 

11.2 

11.2 

Ammonia 

4.3 

7.9 

7.7 

Cyanamide 

4.6 

8.0 

8.0 

Formamide 

4.9 

8.8 

8.8 

Urea 

4.8 

8.0 

7.9 

Ethyl  carbamate 

5.6 

9.0 

8.8 

Acetic  acid 

5.1 

9.1 

9.3 

Oxalic  acid  (o) 

3.2 

6.7 

6.9 

Oxalic  acid  (/ 3 ) 

3.5 

7.3 

7.5 

Succinic  acid 

5.2 

9.1 

9.1 

Adamant  ane 

4.8 

7.5 

7.6 

Hexamine 

5.0 

7.5 

7.8 

Trioxane 

5.9 

9.7 

9.5 

1 ,4-Cyclohexane-di-one 

3.5 

7.0 

7.0 

Benzene 

4.3 

6.8 

6.9 

Naphthalene 

3.2 

5.2 

5.4 

Anthracene 

2.1 

3.9 

4.1 

Pyrazine 

2.8 

6.2 

6.0 

Triazine 

3.0 

6.2 

6.3 

Pyrazole 

4.8 

7.6 

8.1 

Imidazole 

4.8 

7.6 

7.9 

Uracil 

3.4 

6.4 

6.4 

Cytosine 

3.4 

6.1 

6.1 

V 

V 


FIG.  2:  (Color  online)  Absolute  deviations  and  mean  absolute 
deviation  (MAD)  in  the  calculated  quasi-particle  gap  between 
the  OT-SRSH  and  GoWo  for  the  X23  set  of  molecular  solids. 
For  OT-SRSH  calculations,  the  Go  Wo  computed  RPA  macro¬ 
scopic  dielectric  constant  (eRPA)  is  used. 

optical  transitions,  we  only  consider  singlet  excitations. 
The  complete  set  of  optical  absorption  spectra  is  shown 
in  Figs.  3  and  4,  for  solids  based  on  aliphatic  and  aro¬ 
matic  molecules,  respectively.  In  addition,  the  energy  of 
the  lowest  singlet  excitation  and  the  position  of  lowest- 
lying  main  optical  absorption  peak  (i.e.,  the  lowest  sin¬ 
glet  excitation  with  non-negligible  oscillator  strength)  are 
provided  in  Table  II,  along  with  experimental  values  for 
comparison,  where  available. 

It  is  readily  observed  from  Figs.  3  and  4  that  absorp¬ 
tion  spectra  computed  with  the  TD-OT-SRSH  approach 
do  indeed  agree  well  with  those  computed  using  GoWo- 
BSE,  across  the  board,  over  a  range  of  several  eV.  Specif¬ 
ically,  the  position  of  intense  peaks,  found  by  the  two 
approaches,  typically  agrees  within  ~0.2-0.3  eV,  which  is 
as  good  as  can  possibly  be  expected  given  an  accuracy  of 
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Ammonia 


Energy  (eV) 
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Cyan  amide 
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FIG.  3:  (Color  online)  The  imaginary  part  of  the  dielectric  function  of  aliphatic-molecule-based  solids  in  the  X23  data  set, 
calculated  using  Go  Wo /BSE  (black,  solid  lines)  and  TD-OT-SRSH  (red,  dashed  lines). 
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Benzene 


Naphthalene 


Anthracene 


Pyrazine 


Imidazole 


Uracil 


Cytosine 


FIG.  4:  (Color  online)  The  imaginary  part  of  the  dielectric  function  of  aromatic-molecule-based  solids  molecules  in  the  X23 
data  set,  calculated  using  GoVFo/BSE  (black,  solid  lines)  and  TD-OT-SRSH  (red,  dashed  lines). 


~0.1  eV  at  best  for  either  approach  separately.  This  ob¬ 
servation  is  quantified  in  Fig.  5,  which  shows  deviations 
between  peak  positions  using  TD-OT-SRSH  and  GW- 
BSE,  for  the  lowest-energy  peak  (top)  and  for  all  peaks 
shown  in  Figs.  3  and  4  (bottom).  Clearly,  the  mean  abso¬ 
lute  deviation  is  only  0.2  eV  for  either  the  lowest-energy 
peak  or  all  shown  peaks,  with  deviations  rarely  exceeding 
0.3  eV. 

As  expected  based  on  known  shortcomings  of  TDLDA 
[13,  14,  27],  the  TD-OT-SRSH  data  offer  an  improvement 
over  TDLDA  data  that  is  not  only  quantitative  but  also 
qualitative.  This  is  demonstrated  in  Fig.  6  for  two  repre¬ 
sentative  cases  -  the  ammonia  and  uracil  solids.  For  am¬ 
monia,  TDLDA  produces  an  extended  spurious  absorp¬ 
tion  “tail”,  starting  at  ~4.5  eV,  whereas  both  TD-OT- 
SRSH  and  GW-BSE  predict  a  sharp  onset  of  absorption, 
at  ~7.0  eV.  The  same  phenomenon  has  been  previously 
observed  for  a  non-molecular  solid  -  LiF  [29].  For  uracil, 


TDLDA  produces  a  spurious  peak  at  ~3.8  eV.  Further 
analysis  (not  shown  for  brevity)  reveals  that  this  peak  re¬ 
sults  from  LDA  mis-ordering  of  the  HOMO  and  HOMO-1 
orbitals,  which  is  remedied  by  the  OT-SRSH  calculation 
and  removes  the  false  peak. 

In  many  of  the  molecular  solids,  the  lowest  singlet  exci¬ 
tation  possesses  a  small  matrix  element  and  contributes 
little  to  the  optical  spectrum.  This  is  reflected  in  Ta¬ 
ble  II,  where  the  energy  of  the  first-excited  state  is  often 
predicted  to  be  quite  different  from  the  energy  of  the 
lowest-lying  absorption  peak,  using  either  TD-OT-SRSH 
or  GW-BSE.  We  found  that  larger  differences  between 
the  two  methods  often,  but  not  always,  arise  for  these 
low-absorption  excitations,  despite  the  excellent  agree¬ 
ment  in  predictions  of  the  fundamental  gap  and  the  high- 
absorption  excitations.  As  an  example,  for  ammonia  the 
lowest-energy  singlet  excitation  is  predicted  to  be  6.7  or 
6.6  eV  using  TD-OT-SRSH  or  GW-BSE,  respectively, 
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FIG.  5:  (Color  online)  Absolute  deviations  between  the  TD- 
OT-SRSH  and  GoWo-BSE  computed  peak  positions  in  the 
optical  spectra  of  the  X23  molecular  solid  set.  Top:  lowest 
peak  position.  Bottom:  Mean  absolute  deviations  of  all  peaks 
shown  in  Fig.  3  or  4. 

with  both  values  in  excellent  agreement  with  the  experi¬ 
mental  values  of  6.6  eV.  [72]  A  similar  picture  emerges  for 
adamantane.  But  for  urea  or  1,4-cyclohexane-di-one,  the 
difference  between  the  two  predictions  is  a  much  larger 
and  entirely  non-negligible  0.6  eV.  We  note  that  these 
lowest-lying  transitions  often  involve  transitions  between 
highly  localized  orbitals,  which  can  exhibit  large  self¬ 
interaction  errors.  Therefore  the  discrepancy  may  be  due 
to  remaining  issues  in  the  TDDFT  calculation,  but  may 
well  be  also  due  to  LDA  being  an  insufficient  starting 
point  for  the  GW-BSE  calculation.  [82-84]  We  note  that 
such  issues  are  not  necessarily  related  to  strongly  local¬ 
ized  states,  as  both  GW  starting  point  and  remaining 
TD-OT-SRSH  inaccuracies  have  also  been  pointed  for 
benzene  and  oligoacene  molecules  and  solids  [29,  42,  84- 
86] .  As  this  barely  affects  the  optical  spectrum,  the  mat¬ 
ter  is  not  pursued  further  here. 

C.  Effect  of  the  dielectric  constant 

To  facilitate  comparison  to  GW-BSE,  which  relies  on 
evaluation  of  the  dielectric  function  using  the  random- 
phase  approximation  (RPA),  all  OT-SRSH  and  TD-OT- 
SRSH  results  reported  above  were  obtained  using  eRPA. 
Here,  we  explore  the  effect  of  basing  the  calculation  on 
an  evaluation  of  the  dielectric  constant  using  the  in¬ 
expensive  many-body  dispersion  (MBD)  method  [59]. 
A  comparison  of  dielectric  constants  and  quasi-particle 
gaps  obtained  from  using  eRPA  and  eMBD  (not  shown  for 


brevity)  reveals  that  eMBD  >  eRPA  throughout  the  X23 
set,  likely  because  it  is  computed  based  on  PBE,  which 
tends  to  overestimate  polarizabilities.  Therefore,  the 
MBD-computed  fundamental  gaps  are  generally  smaller 
than  RPA-computed  ones.  While  the  difference  is  often 
small,  it  can  be  substantial  -  as  much  as  0.8  eV  for  suc¬ 
cinic  acid  and  0.7  eV  for  uracil.  It  would  be  interesting  to 
examine  whether  an  iteration  of  the  e  value,  from  MBD 
calculations  based  on  the  OT-SRSH  calculation,  would 
result  in  improved  agreement. 

While  the  differences  in  the  dielectric  constant  do  affect 
the  fundamental  gap,  their  effect  on  the  optical  spectra 
is  much  smaller.  This  is  reasonable,  as  the  fundamental 
gap  reflects  charged  excitations,  whereas  optical  excita¬ 
tions  are  neutral.  The  effect  of  the  dielectric  constant  on 
the  TD-OT-SRSH  absorption  spectra  is  demonstrated  in 
Figure  7,  for  the  case  of  the  acene-based  molecular  solids. 
Clearly,  the  effect  of  e  is  marginal  (e.g.,  differences  of  ~0.1 
eV  at  most  in  the  absorption  peak  position  for  the  ben¬ 
zene  solid).  The  effect  on  the  lowest  singlet-excitation 
energy  is  equally  small.  Therefore,  using  MBD  dielectric 
constants  leads  to  an  inexpensive  and  predictive  calcula¬ 
tion  of  optical  spectra  in  molecular  solids. 


IV.  CONCLUSIONS 

In  conclusion,  we  have  computed  fundamental  gaps 
and  optical  spectra  for  the  entire  X23  benchmark  set  of 
molecular  solids  using  the  recently  developed  optimally- 
tuned  screened  range-separated  hybrid  functional  ap¬ 
proach.  In  this  two-stage  approach,  optimal  tuning  of 
a  range-separated  hybrid  functional  is  first  used  for  an 
accurate  and  predictive  calculation  of  the  gas-phase  elec¬ 
tronic  structure.  Dielectric  screening  is  then  built  into 
the  functional  to  obtain  a  self-consistent  prediction  for 
the  solid-state  electronic  structure  and  optical  proper¬ 
ties.  The  obtained  results  have  been  compared  to  many- 
body  perturbation  theory  calculations  within  the  GW- 
BSE  approach.  Agreement  has  been  found  to  be  very 
good  to  excellent  throughout,  with  somewhat  larger  dif¬ 
ferences  possible  for  optically  dark  singlet  excitations 
that  do  not  affect  the  optical  spectrum.  Furthermore, 
we  have  shown  that  inexpensive  evaluation  of  the  dielec¬ 
tric  constant  using  many-body  dispersion  is  sufficient  for 
obtaining  accurate  optical  spectra,  opening  the  door  to 
low-cost,  fully  predictive  calculation  of  optical  spectra  in 
molecular  solids. 

Acknowledgments 

This  work  has  been  performed  in  collaboration  with 
Sivan  Refaely- Abramson  and  Arun  K.  Manna  (Weiz- 
mann  Institute  of  Science),  Anthony  M.  Reilly  (Cam¬ 
bridge  Crystallographic  Data  Centre),  Alexandre  A. 
Tkatchenko  (University  of  Luxembroug),  and  Jeff  B. 
Neaton  (UC  Berkeley). 


DISTRIBUTION  A.  Approved  for  public  release:  distribution  unlimited. 


Ammonia  Uracil 


FIG.  6:  (Color  online)  The  imaginary  part  of  the  dielectric  function  of  the  ammonia  and  uracil  molecular  solids,  calculated 
using  GqWo/BSE  (black,  solid  lines),  TD-OT-SRSH  (red,  dashed  lines),  and  TDLDA  (dotted,  blue  lines). 


Benzene  Naphthalene  Anthracene 


FIG.  7:  (Color  online)  The  imaginary  part  of  the  dielectric  function  for  the  aromatic  acene  molecular  solids  (benzene,  naph¬ 
thalene  and  anthracene),  calculated  using  GoRo/BSE  (black,  solid  lines),  TD-OT-SRSH  based  on  eRPA  (red,  dashed  lines)], 
and  TD-OT-SRSH  based  on  eMBD  (blue,  dash-dotted  lines).  The  different  e  values  are  denoted  in  the  figure. 
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TABLE  II:  Energy  of  the  lowest  singlet  excitation  and  position  of  lowest-lying  main  optical  absorption  peak,  calculated  using 
TD-OT-SRSH  and  the  GoWo-BSE  for  the  X23  set  of  molecular  solids.  Experimental  values,  where  available,  are  provided  for 
comparison.  All  energies  are  in  eV. 


Molecular  Solid 

First  Excited  State  (Si)  Energy 

Optical  Peak  Position 

TD-LDA 

TD-OT-SRSH 

GoWo-BSE 

Expt. 

TD-LDA 

TD-OT-SRSH 

GoWo-BSE 

Carbon  dioxide 

6.4 

8.9 

8.3 

8.9[71] 

7.2 

10.7 

10.8 

Ammonia 

4.3 

6.7 

6.6 

6.6  [72] 

4.6 

7.1 

7.1 

Cyanamide 

4.6 

6.0 

5.4 

- 

5.6 

7.3 

7.1 

Formamide 

4.5 

5.9 

5.4 

- 

5.3 

7.5 

7.5 

Urea 

4.8 

7.1 

6.5 

6. 2  [73] 

6.4 

7.4 

7.1 

Ethyl  carbamate 

5.6 

7.3 

6.5 

- 

7.2 

7.7 

7.5 

Acetic  acid 

5.1 

5.8 

5.6 

- 

5.5 

8.2 

8.3 

Oxalic  acid  (a) 

3.2 

4.7 

4.4 

- 

3.7 

6.0 

6.3 

Oxalic  acid  (yd) 

3.4 

4.8 

4.5 

- 

3.7 

6.1 

6.3 

Succinic  acid 

5.0 

6.2 

5.7 

- 

5.4 

8.7 

8.6 

Adamantane 

4.8 

6.8 

6.8 

6.5[74] 

6.8 

8.2 

8.5 

Hexamine 

4.7 

6.0 

6.2 

- 

5.4 

6.6 

6.9 

Trioxane 

5.8 

8.2 

7.9 

- 

6.3 

8.3 

8.1 

1 ,4-Cyclohexane-di-one 

3.2 

4.5 

3.9 

- 

3.9 

6.2 

6.3 

Benzene 

4.3 

5.4 

4.9 

4.7[75] 

4.9 

6.3 

6.5 

Naphthalene 

3.1 

4.2 

4.1 

3.9  [76] 

3.3 

4.7 

4.4 

Anthracene 

2.0 

2.9 

3.2 

3.1  [77] 

2.1 

3.0 

3.2 

Pyrazine 

2.7 

4.0 

3.5 

3.8[78] 

3.0 

4.0 

3.5 

Triazine 

2.8 

4.4 

3.9 

3.7[79] 

3.4 

4.6 

4.2 

Pyrazole 

4.8 

6.7 

6.5 

- 

5.5 

6.8 

6.5 

Imidazole 

4.5 

6.3 

6.3 

- 

5.0 

6.4 

6.3 

Uracil 

3.4 

4.6 

4.9 

4.5[80] 

3.8 

5.1 

4.9 

Cytosine 

3.4 

4.7 

4.9 

4.4[81] 

3.8 

5.0 

4.9 
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